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ABSTRACT : In this paper, we apply the generalized Bernoulli sub-ODE method to seek exact solutions for 
nonlinear evolution equations arsing in the theory of mathematical physics. For testing the validity of this 
method, we use it to derive exact solutions for the nonlinear Schrodinger (NLS) equation. As a result, some exact 
solutions for it are successfully found with the aid of the mathematical software Maple. 
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I. INTRODUCTION 

Research on solutions of NLEEs is a hot topic in the literature. So, the powerful and efficient methods 
to find analytical solutions and numerical solutions of nonlinear equations have drawn a lot of interest by a 
diverse group of scientists. Many efficient methods have been presented so far. 

Recently, seeking the exact solutions of nonlinear equations has getting more and more popular. Many 
approaches have been presented so far such as the homogeneous balance method [1,2], the hyperbolic tangent 
expansion method [3,4], the trial function method [5], the tanh-method [6-8], the nonlinear transform method 
[9], the inverse scattering transform [10], the Backlund transform [11,12], the Hirotas bilinear method [13,14], 
the generalized Riccati equation [15,16], the Jacobi elliptic function expansion [17,18], the complex hyperbolic 
function method [19-21], the generalized Bernoulli sub-ODE method [22,23] and so on. 

In this paper, we present an application for the generalized Bernoulli sub-ODE method, and seek exact 
solutions for those nonlinear evolution equations arsing in the theory of mathematical physics. 
The rest of the paper is organized as follows. In Section 2, we describe the Bernoulli sub-ODE method for 
finding travelling wave solutions of nonlinear evolution equations, and give the main steps of the method. In the 
subsequent section, we will apply the method to find exact travelling wave solutions for the nonlinear 
Schrodinger equation. In the last section, some conclusions are presented. 

II. DESCRIPTION OF THE BERNOULLI SUB-ODE METHOD 

In this section we present the solutions of the following ODE: 

G '+ XG = /jG 2 (1) 

where X * 0,G = G(i) . 

When jj. * 0 , Eq. (1) is the type of Bernoulli equation, and we can obtain the solution as 

l (2) 

G = V ' 

— + de 
X 

where d is an arbitrary constant. 

Suppose that a nonlinear equation, say in two or three independent variables x , y and t , is given by 

P (u ,U It ,U ,11 ,11 ,u ,u ,u ) = 0 (3) 

where u = u(x, y ,t) is an unknown function, P is a polynomial in u = u(x,y,t) and its various partial 
derivatives, in which the highest order derivatives and nonlinear terms are involved. 
Step 1 .We suppose that 

u(x, y,t) = »(#),£ = f (x, y,t) (4) 

The travelling wave variable (4) permits us reducing Eq. (3) to an ODE for u = u{%) 

P(u,u',u", ) = 0 (5) 
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Step 2. Suppose that the solution of (5) can be expressed by a polynomial in G as follows: 

k(£) = a m G m + a m _ l G m ~ 1 + (6) 

where c = G(£) satisfies Eq. (1), and a m ,a m are constants to be determined later, a m ^ 0 . The positive 

integer m can be determined by considering the homogeneous balance between the highest order derivatives and 
non-linear terms appearing in (5). 

Step 3. Substituting (6) into (5) and using (1), collecting all terms with the same order of G together, the 
left-hand side of Eq. (5) is converted into another polynomial in G . Equating each coefficient of this 
polynomial to zero, yields a set of algebraic equations for a m , a m , . . .X , ju . 

Step 4. Solving the algebraic equations system in Step 3, and by using the solutions of Eq. (1), we can 
construct the travelling wave solutions of the nonlinear evolution equation (5). 

In the subsequent sections we will illustrate the validity of the proposed method by applying it to solve several 
nonlinear evolution equations. 

In the subsequent sections we will illustrate the proposed method in detail by applying it to Fitzhugh-Nagumo 
equation. 

III. APPLICATION OF THE BERNOULLI SUB-ODE METHOD FOR NLS EQUATION 

In this section, we will consider the following NLS equation: 

^ ; -^„ + 2(|^| 2 -p 2 )^ = 0 (7) 

where <j> is complex wave function and p is a constant. 

Since cj> = </>(x,t) in Eq. (7) is a complex function, we suppose that 

^ = «(|)exp[i(ax + fit)], £ = k (x + 2a t) (8) 
where the constants a , fi ,k can be determined later. 

By using (8), (7) is converted into an ODE 

(-fi + a 2 - 2p 2 )u + 2u - k 2 u " = 0 (9) 
Suppose that the solution of (10) can be expressed by a polynomial in G as follows: 

«(<?) = £ a,G' 



(10) 



wherea. are constants, and G = G(£) satisfies Eq. (2.1). 

Balancing the order of u 3 and u " in Eq. (9), we obtain that 3m = m + 2=> m = 1 .So Eq. (10) can be 
rewritten as 

k(£) = fljG + a 0 ,a 1 * 0 (11) 

a l , a 0 are constants to be determined later. 

Substituting (11) into (9) and collecting all the terms 

with the same power of " together, the left-hand side of Eq.(9) is converted into another polynomial in G . 
Equating each coefficient to zero, yields a set of simul-taneous algebraic equations as follows: 

G° : -/?«„ + a 2 a 0 - 2p 2 a 0 + 2a 0 ' = 0 

1 2 2 2 2 2 

G : - p a 1 - 2 p a x + a a t - k a t X + 6 a a x = 0 
G 2 : 6 a 2 a 0 + 3k 2 Xa { p - 0 
G 3 : 2a' - 2k 2 a % p 2 = 0 

Solving the algebraic equations above, yields: 



Case 1: 



1 (12) 



a, = -k fj ,a 0 = —kX,/3 = a' - 2 p + 2 a 
2 
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Substituting (14) into (13), we have 

= -kpG + -kX > <? = k(x+ 2at) ■ (13) 



2 



Combining with Eq. (2.2) and considering u = v " _1 , we can obtain the travelling wave solutions of (7) as 
follows: 

ii 222 (14) 

= [-kfi( ) + -U]exp[i(ai + (a ~2p + 2a„ )»)] 



^ . , H 2 



X 

where d is an arbitrary constant. Then we have 



11 2,2 (15) 

0 t (x,t) = [-kfi( ) + —kl]exp[i(a x + (a -2p~ + 2a a )t)] 



Case 2: 



■ de 



a, = k/u,a 0 = - — kX, B = a 1 - 2 p 1 + 2a 0 2 
2 



(16) 



Substituting (14) into (13), we have 



u 2 (%) = kfiG --kX <4 = k(x + 2at) ■ (17) 
2 

l 

Combining with Eq. (2.2) and considering u = v " _1 , we can obtain the travelling wave solutions of (7) as 
follows: 

I 1 2 2 2 (18) 

0>,(f) = [kfi( )--kX]exp[i(ax + (a -2p +2a )t)] v ' 

M , h 2 
— + de 

X 

where d is an arbitrary constant. 
Then we have 

II 222 (19) 

0 2 (x,t) = [-k fi( ) + — kX]exp[i(ax + (a -2p + 2a Q )»)] 

M , Ikix+lat) 2 

— + de 
X 

IV. CONCLUSION 

In the present work, we have found some new travelling wave solutions for the nonlinear Schrodinger 
(NLS) equation by the Bernoulli sub-ODE method. As one can see, this method is straight ward, and can be 
fulfilled with the aid of the mathematical software Maple. Being concise and simple, this method is one of the 
most effective approaches handling nonlinear evolution equations, and can be used to many other nonlinear 
problems. 
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